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ABSTRACT 

Recent ideas for the origin and persistence of the warps commonly observed in disc galaxies 
have focused on cosmic infall. We present A^-body simulations of an idealized form of cosmic 
infall onto a disc galaxy and obtain a warp that closely resembles those observed. The inner 
disc tilts remarkably rigidly, indicating strong cohesion due to self-gravity. The line of nodes 
of the warp inside i?26.5 ~ 4.5i?d is straight, while that beyond i?26.5 generally forms a 
loosely-wound, leading spiral in agreement with Briggs's rules. We focus on the mechanism 
of the warp and show that the leading spiral arises from the torques from the misaligned 
inner disc and its associated inner oblate halo. The fact that the line of nodes of most warps 
forms a leading spiral might imply that the disc mass is significant in the centre. If the line 
of nodes can be traced to very large radii in future observations, it may reveal information on 
the mass distribution of the outer halo. The warp is not strongly damped by the halo because 
the precession rate of the inner disc is slow and the inner halo generally remains aligned with 
the inner disc. Thus even after the imposed quadrupolar perturbation is removed, the warp 
persists for a few Gyrs, by which time another infall event can be expected. 
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1 INTRODUCTION 

The optically visible parts of galactic discs are usually remarkably 
thin and flat, whereas the more extended HI discs of many edge- 
on galaxies appear noticeably warped with an i ntegral sign shape. 
Stella r warps do exist in some galaxies (e.g. iReshetnikov et alJ 
|2002|) . but are much less pronounced than the warps in the ex- 
tended gaseous discs. We infer that warps are a gravitational phe- 
nomenon because weak stellar w arps, if present, follow the same 
warped plane as the gaseous ones ICox et alll996l) . 

Warps are extre mely common - at least half of all spiral 
galaxies are warped: Bo smal i 19911) found 12 out of 20 edge-on 
_spiral galaxies are warped. The more recent HI observations by 
lOarcia-Ruiz et alj 120021) found that 20 out of 26 edge-on galax- 
ies are warped. Since these surveys detected only warps whose line 
of nodes is close to ou r line of sight, the true frac tion of warps must 
be even higher. In fact lGarcia-Ruiz et alj f 2002) also found that all 
galaxies possessing an HI disc more extended than the optical one 
are warped. The ubiquity of warps suggests that warps are either 
repeatedly regenerated or long-lived. 

Warps can be detected kinematically even when the system is 
not edge-on. iBriggsH 1 9901) studied a sample of 12 warped galaxies 
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with high-quality 21 -cm data, and found that galactic warps obey 
some fairly simple rules: 

(i) the HI layer typically is coplanar inside radius R25 , the radius 
where the B-band surface brightness is 25 mag arcsec"^ — 25/iB -^ 
and the warp develops between R25 and R26.5 (aka the Holmberg 
radius); 

(ii) the line of nodes (LON) is roughly straight inside -R26.5; 

(iii) the LON takes the form of a loosely-wound leading spiral 
outside R26.5- 

The origin and per sistence of warps still presents a puzzle. 
iHunter & Toomrd il969h . before the discovery of dark matter ha- 
los, studied the bending dynamics of isolated thin discs. They found 
that discrete warp modes in a cold, razor-thin disc do not exist un- 
less the edge is unrealistically sharp. Such modes (standing waves) 
can be realized by the superposition of outgoing and ingoing waves, 
provided that waves incident on the edge can reflect. Ho wever, a 
disc with a smooth edge does not reflect bending waves JToomrd 
1983) . 

Subsequent ideas of warp formation rely in some way 
on the interaction between the disc and its dark matter halo. 
Dekel & Shlosman ( 1983) and Toomi'e ( 1983) suggested that a flat- 
tened halo misaligned with the disc could form a long-lasting warp. 
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If the disc is exponential with scale-length R^ and central surface bright- 
ness /gg ^ 21.7 mag arcsec"^ = 21.7/^b (Freeman 1970, the so-called 
Freeman's law), R25 ^ S.OiJd and -R26.5 ^ AAR^. 
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ISparke & CasertandiigsS) and lKuiikeiJJ199lh invoked a rigid flat- 
tened halo and obtained persisting warps (dubbed "modified tilt 
modes"), which were insensitive to the details of the disc edge. 
[Lovelace ( 1998) studied the tilting dynamics of a set of rings also 
in a rigid halo, but generally assumed that the inner disc lay in the 
symmetry plane of the spheroidal halo. 

Dark matter halos are not rigid, ho wever, and a responsive halo 
alters the dynamics in several ways. iNelson & Tremain e ( 1995) 
showed that the precession of a misaligned disc inside a flattened 
halo would sap energy from the warp mode through dynamical 
friction, damping it on timescales much shorter than a Hubble 
time. Numerical studies using A'^-body halos als o found the discrete 
warp mode does not survive JPubinski & Kuiike n ( 1995) f ound that 
the inn er halo and disc quickly settle into alignment. .Binnev et alj 
il998h confirm that a discrete warp mode in a rigid halo does not 
survive in a simulation with a responsive halo, and conclude that 
the inner halo could never be significantly misaligned from the in- 
ner disc^ 

lOstriker & BinnevI il989|) drew attention to the likely mis- 
alignment with the disc axes of late infalling material in hierar- 
chical galaxy formation models, and proposed that warps arise due 
to the slewing of the galactic pot ential as material with m isaligned 
angular momentum is accreted. lOuinn & BinnevI il992h showed, 
for a realistic cosmic scenario, that the mean spin axis of a galaxy 
must slew as late arriving material rains down on the early disc. The 
less-than-critical matter density in modern ACDM universe models 
implies that infall is less pervasive at later times, but it manifestly 
continues to the present day in gravitationally bound environments 
such as that of the Local Group. 

Jiang & Binnev ( 1993 hereafter JB99) present results of an 
experiment in which a disc is subjected to the torque from a mis- 
aligned, massive torus at a large radius. This well-defined pertur- 
bation is a very crude model of an outer halo that is rotationally 
flattened, with a spin axis misaligned with that of the disc. It is 
misaligned and farther out because, in hierarchical scenarios, the 
mean angular momentum of the later arriving outer halo is proba- 
bly misaligned from that of the original inner halo and disc. They 
concede that the accretion axis is in reality likely to slew continu- 
ously over time, so a model with a constant inclination is somewhat 
unrealistic. 

Here we use high quality A'^-body simulations of a model of 
this type. We improve on the work by JB99 in the following main 
aspects: 



• We use a disc of particles with random motion, whereas JB99 
employed a disc composed of rigid rings coupled only by grav- 
itation. Thus JB99 did not include random motion of the disc 
stars, which has been shown to add to the stiffness of the disc 
JDebattista & SellwooJl999l) : 

• We use a much more extensive disc so that we can study the 
warp - JB99 truncated the disc at 4_Rd. The behaviour of the LON 
beyond this radius is important for comparing our results with ob- 
servations of warps (e.g. Briggs's rules); 

• In il6.2l we remove the forcing perturbation at later times and 
follow the evolution of the distorted disc fully self-consistently. 



We obtain flat inner discs and long-lasting warps in the outer disc 
that match all of Briggs's rules quite well. We present a detailed 
analysis of the dynamics and show, for example, that the persis- 
tence of warps is not nearly as perplexing as is currently believed. 



Figure 1. A sketch of initial model setup. The dark matter halo is not shown. 



2 MODEL SETUP AND SIMULATION DETAILS 

Following JB99, our simulations include three distinct mass com- 
ponents: a disc, a halo, and an accreted torus. The plane of the torus 
is inclined to that of the disc, as sketched in FigureQ 

Our simulations are set up as follows: The initial disc has the 
exponential surface density profile: 



^(^) = T^ exp 



2nRl 

where R is the cylindrical radius and Rd and Md are the scale ra- 
dius and the total mass of the disc, respectively. We truncate the disc 
at i? = 8Rd, spread the particles vertically as the isothermal sheet 
with a locally-defined equilibrium vertical velocity dispersion, and 
embed it in the halo; the rms vertical thickness of the initial disc is 
0.1i?d at any radius. We set the initial velocities of the disc parti- 
cles such that the Toomre stability parameter Q ~ 1.5 and the disc 
is in rotational balance. 

The halo is set up so as to be in equilib rium with the disc using 
the procedure described in Appendix A of iDebattista & Sellwood 
12000.) . The distribution function has a King-model form with 
$(0)/cr^ — 2.0, where a is the central vel ocity dispersion and 
'i'iO) is the relative potential at the centre I'B innev & Tremaina 
[1.987), which in our case includes a contribution from the disc. We 
set the tidal radius rt = 18i?d. and total mass — 9.0Md- The 
King radius is Vc ~ 5.6/Jd and the half-mass radius of the halo is 
rh ^ 5Rd- The inner halo is mildly oblate initially due to the grav- 
ity of the disc. The halo we adopted here is massive, as we attempt 
to use a similar halo profile used by JB99 in which A/h/A/d ~ 5.0 
for r < 3.2Rd. Other experiments with halos having different mass 
profiles, to be reported elsewhere, have revealed that the conclu- 
sions presented here are not strongly affected by the choice of halo; 
warping behaviour is less sensitive to the density profile or total 
mass of the halo than to the masses and sizes of the disc and torus. 

Following JB99, we gradually inject particles into an inclined 
torus, causing its mass to increase linearly from at i = 0, to 
Mt = 2.5Md at t = 200. The torus, which has radius Rt, 
is uniform in azimuth and has Gaussian density profile in cross- 
section, with a standard deviation b; in our case, Rt — 15.0i?d and 
b = 0.2_Rd. The torus plane is inclined at 15° relative to the ini- 
tial disc plane and we adopt the line of nodes as the j/-axis of our 
coordinate system. 

The initial velocities of injected torus particles are equal to the 
local circular speed in magnitude and directed tangentially around 
the torus, as in JB99. Thus the torus has very large angular momen- 
tum, and therefore does not precess significantly from its original 
plane in response to the torque from the disc. However, a cold torus 
is dynamically unstable and breaks into radially extensive spiral fil- 
aments after a few Gyrs when allowed to evolve self-consistently. 
The breakup of the cold torus depends on many factors that we do 
not fully understand — e.g., tori with the same parameters dissolve 
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Table 1. Numerical parameters used in the canonical simulation 



Cylindrical grid Spherical grid 



Grid size {Nii,N^,N^) 

= (102, 128, 125) 

Angular compnts .... ^ m ^ 8 

Outer radius (_Rd) ... 8.0 

2-spacing (i?d) 0.02 

Softening length (R^) 0.02 

N O.IM 

Time step 0.04 



Nr 



:401 



^/ ^4 
18.0 



0.9M 
0.04 



differently and at different times in runs with rigid and live halos. 
In order to avoid the inconvenient complication of time-varying, 
non-axisymmetric perturbing forces acting on the disc, we gener- 
ally keep the torus particles fixed in their initial positions relative 
to the galactic centre. This stratagem enables us to compare the dy- 
namical evolution of warps in different experiments with the same 
torus parameters. 

All 1.25 million particles belonging to massive components in 
our simulations have equal masses. 

We adopt R^ and M^ as our units of length and mass, 
respectively, and our time units are therefore dynamical times 
(R^/GMd)^^^. Generally all quantities are expressed in units such 
that G = A/d — Rd ~ 1 unless otherwise noted. These units can 
be scaled to physical values as desired; we adopt one possible scal- 
ing, choosing R^ = 2.5 kpc and a unit of time of 10 Myr implies 
Md = 3.47 X 10^° Mq; our unit of velocity scales to 244 km s~\ 

W e use the hybrid particle-mesh scheme described in de- 
tail in ISellwood i2003l Appendix B). The self-gravity of the 
disc is computed on a high-resolution cylindrical polar grid (see 
[Sellwood & Valluri 1997; Shen & Sellwood 2004), while that of 
the halo is computed using a surface harmonic expansion on a 
spherical grid. We have carefully checked that the preferred plane 
of the cylindrical polar grid does not introduce artificial forces be- 
tween a tilted disc and the grid, which requires a fine polar grid. 
Our adopted numerical parameters are summarized in TableQ 

Over a period of 400 dynamical times, energy is conserved to 
about 0.02%, angular momentum components (Lx, Ly, Lz) change 
no more than 1%, and the absolute change of linear momenta (px, 
Py, pz) are less than 0.06 in units we adopted. Without the mis- 
aligned torus, these global integrals are conserved with a precision 
almost one order of magnitude better. 

Our main warp diagnostic is the tip-LON diagram introduced 
by|Bri2£s ( 1990). We construct it from an analysis of the disc par- 
ticles, which are binned into annuli. The inner disc is very stiff 
and remains very closely coplanar within the innermost 3Rd', we 
therefore use a single bin for all particles inside this radius and 
adopt cylindrical bins (spherical bins are preferred if run for a 
very long time) of equal radial widths 0.7-Rd outside this radius 
lo R — 7.9/Jd- (The bins remain aligned with the original disc 
axes.) We then compute the eigenvectors of the inertia tensor to 
determine the orientation of the disc element, 9 and 0lon in each 
bin. We use the orientation of the inner bin to define the plane of 
the inner disc. We have experimented with various bin widths and 
found this binning scheme gives us the best compromise between 
spatial resolution and noise in estimation of the disc orientation (es- 
pecially for thicker/"hotter" inner discs). We describe the diagram 
and its meaning in il3.3l 
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Figure 2. (a) The warp at i = 400 in our simulation; its morphology closely 
resembles the observed HI warp of NGC 4013 (Bottema 1996, reproduced 
with permission), shown in (b). The length unit shown is the scale length 
Rd of the exponential disc. Note that we have oriented the model so that 
the inner (R < 3Rd) disc lies in the x-y plane, which is perpendicular to 
the paper. 





R/Rd 

Figure 3. (a) Disc tip angle as a function of radius at intervals of 40 time 
units (t = to t = 320 from bottom up), (b) as for (a), but the warp angle 
9 relative to the inner disc. The inner disc is defined by the radial range 
< R < 3Rd. 
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Figure 4. Tip-LON diagrams showing the time evolution of the warp in our 
canonical experiment. The various points indicate 8 (radial coordinate) and 
4>LON (azimuthal coordinate) averaged over annular bins having widths of 
0.7i?d- The warp angle, 6, is larger for annuli with larger radii; the outer- 
most point is for the annulus centered on ij = 7.55R;i- The large circles 
mark warp angles of 10° (dotted) and 20° (solid). 



3 A CANONICAL SIMULATION 

Without the torus, the disc in our model heats mildly; Q rises from 
~ 1.5 to ~ 1.8 over 300 dynamical times, the disc thickens by 
~ 20% and, of course, shows no signs of warping. The virial ratio 
stays at 0.5, confirming the equilibrium of our initial model. The 
disc remains unbarred until very late times. 

We grow the inclined torus in our canonical experiment, in 
which the initially thin and flat disc experiences a growing torque. 
Since the disc is spinning rapidly, it tends to precess about the sym- 
metry axis of the distant torus; differential precession will make the 
disc warp. We present a quantitative discussion of the precession in 
S|4| The much larger moment of inertia of the accreted torus implies 
that its precession in response to the torque from the disc is mini- 
mal, and we therefore keep the torus particles frozen in their initial 
positions relative to the galaxy centre in order to provide steady 
forcing. Some differences in the disc behaviour between this case 
and one in which the torus particles move begin to appear after 
t ~ 320. 



perfect; e.g., our simulated warp starts at around 5Rd whereas the 
warp of N GC 4013 starts t o develop at the edge of its optical disc, 
about 4i?d I DOttemal99D). We could adjust the density profile of 
the halo, which affects the radius at which the disc starts to warp, 
but have not attempted to tailor our model to improve the match. 

Note also that the particles in our simulation are collision- 
less whereas the warp is observed in the dissipative gas compo- 
nent. Nevertheless, the comparison should be reasonable because 
the precession rate of the simulated warp is low - collisional dissi- 
pation in a gas layer be comes impor t ant only when the d ifferential 
precession rate is high jGunnlll979tlTohline et alJll982h . Further- 
more gas often condenses into very compact clouds, which behave 
quasi -ballistically in a similar manner to a stellar system JBinnevI 



3.2 Warp angle 

Figure|3ja) shows &{R), the tip angle relative to the z-axis of the 
initial disc, at various times. In Figure Isjb) we transform to warp 
angles 6, which are measured relative to the orientation of inner 
disc (R < 3Rd) at each time. The warp amplitude at later times 
in Figurelsjb) seems to compa re well with tha t in NGC 4013, both 
visually and quantitatively ("cf . Ib ottemal 199a. Figure 6). 



3.3 LON of the warp Is a leading spiral 

Consistent with our convention for warp angles, we use (Jj-lon for 
the azimuth of the line of nodes between the plane of the inner disc 
and that of the annular bin, and "I"lon for the corresponding line of 
nodes with the initial unperturbed disc plane. 

Figure|4lsho ws the evolutio n of LON of the warp, in the form 
of tip-LON plots iBrigg3ll990l) . The points (small open circles) 
indicate the angle (6, cf>hON} of the symmetry axis of an annulus 
of disc material, relative to the inner disc plane at any given time. 
The radial coordinate is the warp angle (inclination), 0, and the 
azimuthal coordinate indicates the azimuth, ^lon, of the best-fit 
plane to a disc annulus. Points at successively larger ring radii are 
joined sequentially by straight line segments. Note that since the 
warp angle, 6, increases with radius (Figure |3j, the plotted points 
refer to radially ordered annuli from the centre out. 

Figure |4| shows that the LON of the warp precesses in the ret- 
rograde direction, while curving gradually into a leading spiral (i.e., 
the LON advances in the direction of galaxy rotation for succes- 
sively larger radii), consistent with the third rule of Briags ( 199(|). 
We have run many other experiments with various parameters and 
found that the LON to 8Rd always forms a leading spiral when the 
disc is perturbed by an outer inclined torus. The reason for the lead- 
ing spiral is explained in detail in the next section (i|4j. Also the 
LON is fairly loosely wound throughout our simulation. Ib ottemal 
lll99d) concluded that the LON of NGC 4013 winds modestly, by 
about 20°, which is consistent with that in our simulated warp, in 
addition to the nice visual resemblance in Figure|2| 



3.1 Morphology of warps 

Figure |2| shows the edge-on projection of the warped disc in the 
canonical simulation at i = 400, or 200 dynamical times after torus 
growth ended. To reveal the warp more clearly, we have rotated the 
particle coordinates to the frame of the inner disc. 

To illustrate that the warp in our simulation indeed resembles 
those observed, we c ompare it to th e strongly warped hydrogen 
layer of NGC 4013 JB ottemal 1 19961) . The correspondence is not 



4 IN-DEPTH ANALYSIS ON THE SHAPE OF LON 

Since there are many factors that affect the warp morphology, we 
break the dynamics into parts. We first study the simplest aspect an- 
alytically, and then gradually add in separate pieces of the physics 
to illustrate their effects. These simplified experiments are very 
helpful for developing understanding of warp dynamics and the 
reason for the curved LON. We stress that our adopted massive 
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Figure 5. The retrograde precession of the angular momentum of a disc 
annulus under the torque of the accreting torus. The precession is about the 
symmetry axis of the torus. 



torus at a very large radius has so much angular momentum that its 
precession due to the torques from the disc and flattened inner halo 
can be neglected. 

4.1 A test particle disc In a rigid halo 

We first study a very simple model, in which the halo is both rigid 
and spherical. We also replace the self-gravitating disc with mass- 
less test particles and represent only the torus with massive par- 
ticles. This model is sufficiently simple that the behaviour can be 
predicted analytically and provides a welcome check of our numer- 
ical code. We expand the potential of the distant torus as a series of 
multipoles and calculate the disc precession rate analytically, both 
to second- and to fourth-order, in Appendix]^ we give exact for- 
mulae for B and ^lon in Appendix IbI 

We use 20 000 test particles to represent the massless disc and 
adopt a rigid, spherical King-model halo with \['(0)/cr^ = 6.0, 
ft — 18.0, and total mass 9.0. Only the torus is composed of mas- 
sive particles, and is grown from Aft = to 2.5 within 200 time 
units, as before. The inclination angle of the disc with respect to the 
plane of the accreting torus is i = 15° . 

The disc is composed of many concentric rings of massless 
particles, which precess at fixed inclination in a retrograde fashion 
under the torque of the distant accreting torus. The precession is 
illustrated in Figure|5| 

We use the following notation: u^ is the nodal precession rate 
of a narrow annulus of the disc at radius R\ O is the angle between 
the symmetry axis of this annulus and the z-axis (i.e., the normal 
to the initial unperturbed disc); "1?lon is the azimuthal angle of the 
symmetry axis of the annulus projected in the x-y plane, which is 
the plane of the initial unperturbed disc. 

Each annulus precesses through an angle 



A* 



P - / ^p{t) dt, 



(2) 



where the precession rate due to the quadrupole term only (Equa- 
tion lA9t is 

3 GMtR 



^p, torus 



4 R^Vc 



■ cosi (X R for Vc ~ const. 



Clearly A$p increases with R. 

From Equations <B2t and <B5> . we find 

e « 2i sin(A<l>p/2) for small 9, i 

and 



(3) 



(4) 




Figure 6. The angles "1'lon (Irft) and (right) measured from the test par- 
ticle disc discussed in i|4. II Particles at r = 3.15 are shown in the top row 
and at r = 7.05 in the bottom row. The light solid lines are the angles mea- 
sured from 8 particles at each radius; the heavy solid lines are predictions 
based on the fourth-order approximation of a;p (Eguation lASI : the dashed 
lines are results based on the second-order approximation (Eauation lA9t . 
The fourth-order approximations are almost invisible because they overlap 
the curves measured from the 8 particles, indicating excellent agreement, 
while the second-order approximations are clearly inadequate at large radii. 
Note that ^lon has been set to —90° for < 1°, to avoid large scatter 
at early times. 



<E>LON ~ — 7r/2 — A<i>p/2 for small i 



(5) 



For test particles at a given radius, the tip angle O and "1>lon 
can be measured from the orientation of their angular momentum 
vectors. The angles, "I>lon and B, measured from particles at two 
distinct radii are shown in Figure |6l each particle in an annulus 
has its own pair of angles at each instant, which vary systemati- 
cally with phase giving rise to the spread. The mean is in qual- 
itative agreement with these second-order predictions, shown by 
the dashed curves while the fourth order (Eauation lA8> predictions 
(solid curves) agree very well. Since ^lon can have any value 
when O ~ 0, we set ^lon to —90° for O < 1°, to suppress 
noise in these figures. 

The increasing rate of retrograde precession with radius 
causes the warp angle 6 and 0lon to increase outwards also. We 
observe this behaviour in our test-particle simulation as shown in 
Figure0 In this simplified model, the warp appears as the result of 
the differential precession and the LON forms a trailing spiral. 

This simulation also serves as a nice check of parts of our 
numerical code: forces from the torus particles and the time inte- 
gration of the test particle disc are calculated accurately. 



4.2 Including disc self-gravity 

We separate a self-gravitating disc into an inner disc and an 
outer disc. The dense inner disc, which contains most of the 
disc mass, tilts rigidly as a whole; it is strongly coupled together 
due both to self-gravity and to radial velocity spread of stars, 
which communicate stress ac ross the disc via epicyclic excursions 
JDebattista & Sell wooc^l999^. The gradual locking of the disc due 
to self-gravity was also found in Lovelace ( 1998). The low-density 
outer disc, on the other hand, is much less cohesive because self- 
gravity is weak and epicycles are small; the outer disc behaves more 
as a collection of test particles therefore. 

As the torus mass rises, the whole disc starts to precess as de- 
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Figure 7. Tip-LON plots of the test particle disc model in different refer- 
ence frames, (a) The top two rows show © and 'I'lon with respect to the 
fixed coordinate system of Figure|^ (b) The bottom two rows show the an- 
gles 8 and ipLON relative to the inner disc frame at each time, which is more 
appropriate for comparisons with observations. The radial binning scheme 
is the same as that used in Figure |4] but the concentric dotted circles are 
drawn at intervals of 10° . 



scribed in i|4. II The torque from the torus causes the inner disc, 
which is strongly cohesive, to precess slowly as a whole in a retro- 
grade manner about the symmetry axis of the torus, while the outer 
disc precesses more rapidly. The developing misalignment between 
the inner and outer disc causes the particles in the outer disc to feel 
an additional torque from the massive inner disc. The rate of preces- 
sion of the outer disc due to the torque from the inner disc follows 
a different rule from that due to the torus. 

The disc potential at larger radii can be approximated as the 
sum of a monopole and a quadrupole term (.Binnev & Tremaine. 
Il987l Equation 6-84), 



-I'd 



GA4i 



G 

2r3 



P2(c 



2nR'^T.{R')dR' . 



Using equation Q for an exponential disc, we have 



<E>d 



GA/d GMd 37?^ 



P2(cos( 



(6) 



(7) 



In the same manner as for the torus (Appendix |^ see also 
[Kahn & Woltier 1959 ; Gunn 1979) , we find the precession rate of 
test rings due to the disc is 




Figure 8. Each row shows part of the warp evolution for different disc 
masses: (a) M^ = 0.1, (b) A/j = 0.3, (c) Afj = 1.0, and (d) Afj = 2.0. 
Note the A/j = case is shown in Figure0b). The spiral formed by the 
inner LON clearly changes from trailing to leading as A/j is increased. 
Dotted circles mark 10° intervals. 



^p, disc ' 



GA/d 



irVc) 



While still producing retrograde precession, the rate from this 
additional forcing term varies with radius in the opposite sense of 
that from the torus. In the case of forcing by the inner disc, the pre- 



cession rate Equation Js} decreases with increasing radius, which 
in isolation would tend to make the LON develop a leading spiral. 

The outer disc is now subject to two torques, from the torus 
and from the inner disc, both of which cause retrograde precession, 
but at rates that vary radially in opposite senses. Whether the re- 
sulting spiral in the warp LON should lead or trail over the range of 
interest depends on relative magnitudes of the torques, and is best 
determined from simulations. 

We have therefore performed a series of simulations with in- 
creasingly massive discs to show that the torque from the inner disc 
dominates in our canonical simulation, which develops a leading 
spiral in the LON (see Figure|4}. Figure|8|presents the tip-LON di- 
agrams from four experiments with increasingly larger disc masses 
(from A/d = 0.1 to 2.0) to demonstrate the effect of the disc mass. 
These experiments continue to employ the same spherical rigid 
halo and torus parameters as above. The general trend confirms that 
the leading spiral gradually dominates the inner LON as A/d is in- 
creased; the LON has become completely leading in Figure |8jd) 
where A/d ~ 2.0. (Note that the result for Afd = 1 in a rigid halo 
differs from that in Figure|4]where the halo is responsive.) 

4.3 A test particle disc in a live halo 

for Vc ~ const. (8) The rigid spherical halo in the above experiments, cannot respond 
to the precessing disc, and therefore does not affect how the LON 
curves. A halo composed of live particles should acquire a distorted 
shape in response to the fields of the torus and the disc. These ef- 
fects are both included in a fully self-consistent simulation. Here 
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Figure 9. The evolution of the warp for a test particle disc in a live halo. The 
only difference from FigureQis that the live halo becomes slightly oblate 
under the influence of the distant torus. 




Figure 11. The spirality of LON twists from leading to trailing at very large 
radius (the transition occurs near the eighth point from the centre, i.e., the 
annulus centered on R = 7.55Rd)- These tip-LON plots are almost same 
as those at t = 200 and 400 in Figure|4] except that 14 radial bins are used 
to cover all particles, including the additional massless ones; the binning 
scheme is similar to that described in i|3.3l Note the disc rotates in a counter- 
clockwise sense. 



$LON ~ const 
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Figure 10. Schematic illustration showing two possible spiralities of the 
warp LON. (a) A leading spiral caused by a precession rate ujp that is a 
monotonically decreasing function of R such as a;p ~ (.iJp, disc oc iJ"''; 
(b) A trailing spiral resulting from a precession rate uip that is a monoton- 
ically increasing function of R such as uip ^ uip, torus oc R.. The arrows 
outside the plots indicate the direction of galaxy rotation. 



we study these effects separately and demonstrate that the live halo 
does not alter the conclusions we reached in i|4.1l and [4.2l 

In a further experiment, we replace the rigid halo in the model 
of ii4. II with a live one with nearly the same potential. Figure |9| 
shows that in this experiment the LON becomes trailing, as we 
found in i|4. II where the halo is spherically rigid. This again sup- 
ports our argument that the massive disc, instead of other factors 
like a live halo, is the primary reason why the LON forms a leading 
spiral for warps induced by cosmic infall. The small differences in 
the angle of precession between Figures |9| and 0b) indicate that 
the oblateness of the halo (induced by the torus) contributes only 
slightly to the precession of the disc. 

We discuss the response of a live halo with a massive disc in 



m 



4.4 Spirality twist of LON at large radius 

In Figure[8lc), we can clearly see the spirality of LON twists from 
leading to trailing in the tip-LON plot at i = 500. Here we show 
that such spirality twist is a generic feature for warps caused by an 
external quadrupole field giving outward-increasing rate of retro- 
grade precession. 

As shown above, the sense of spirality of the LON is deter- 
mined by the radial variation of the precession rate cup, and is mea- 
sured by $LON in the tip-LON plots. From Equation |2| and |5] we 
have 



- / LUp{R, t) dt 



(9) 



the negative sign indicates that the precession is retrograde. 
Whether the spiral leads or trails depends on whether tjp falls or 
rises as R increases, as illustrated in Fieure fTol 

We have argued that the torque due to the massive inner disc, 
which causes Wp ~ ojp, disc c< i?^* to dominate, causing the 
leading spiral found in our simulations. Beyond some critical ra- 
dius, 7?tw, precession is dictated by the accreting torus, u)p ~ 
<^p, torus oc R, and a trailing LON should be expected. The crit- 
ical radius i?tw can be estimated as follows: The magnitude ratio 
of the torques exerted on a test ring at radius R, due to the torus 
(the first term of Equation I A7I as the second-order approximation) 
and to the massive inner disc (Equation|8}, is approximately 



Tto 



Tdii 



3 GKh „2 

4 i?? 



9 GMd Rl 
2 R R^ 



1 Mt R^ 
6 Afd RlR^ 



The radius of the spirality twist is approximately where 



Tto 



Tdk 



-R=Ht» 



So we have 



Rt' 



6 



mr-i-^"^ 



311/5 



.(10) 



(11) 



(12) 



For our torus parameters (Mt = 2.5, Rt — 15), we find i?tw ~ 
6.1-Rd- In Figure|8jc) where A/d ~ 1-0, the spirality twist occurs 
near the sixth point from the centre (i.e., near the aimulus centered 
on i? = 6.15-Rd), so it agrees well with the expected 7?tw 

i?tw can also be understood in t he context of a Laplacian 
surface, where the net torque vanishes JBinnev & Tremaindl 19871 
chap. 6.6). It is easy to show that Rt„ is also the transition radius 
of the Laplacian surface; for R <C i?tw, the Laplacian surface al- 
most coincides with the inner disc, which dominates the precession 
for small radii, while for R S> -Rtw, the Laplacian surface follows 
the equatorial plane of the torus, which dominates the precession at 
large radii. 

The inner part of a responsive halo becomes oblate due to the 
massive inner disc. The co-aligned oblate inner halo adds signifi- 
cantly to the torque from the inner disc, causing i?tw to be larger 
than estimated above, to perhaps 77?d ^ ^tw ^ 8Rd- So the LON 
spirality twist could occur at radii > 2Rho- 

Thus the spirality twist can be seen only if there is a tracer at 
very large radii. We do not see the spirality twist in Figure|4]because 
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very few particles lie beyond ~ SRd- To reveal the expected spi- 
rality twist, we add into our canonical simulation 20,000 massless 
test particles in the radial range of 3-Rd < R < 12Rd, to mimic 
the outer gas layer. We used 14 radial bins to cover all particles; the 
binning scheme is similar to that described in il3.3l The spirality 
twist shows up in Fi sure fTTI near the seventh from outermost point 
(i.e., near the annulus centered on _R = 7.55-Rd). As expected, the 
spirality twist grows more pronounced with time. 

Note that a twist can occur only if the forcing causes an 
outward-increasing rate of retrograde precession, as for our sim- 
plified torus. More realistic perturbations may not have this effect. 



5 THE HALO RESPONSE 

In this section, we return to our canonical simulation (jsj with the 
live halo and the massive disc (il/d = 1.0). 



5.1 Disc-Inner halo alignment 

Figure fT2l shows the orientation of the halo in three spherical bins 
(0 < r < 4i?d, 4i?d < r < 8i?d, and 8Rd < r < 127?d) 
measured from the canonical simulation. The iimer halo (heavy 
line) follows the inner disc (dashed) closely, while the outer halo 
(light line) quickly aligns with the outer torus, which has a tip an- 
gle Q ~ 15° and azimuth 3>lon ~ 180°. The halo orientation in 
the middle bin is intermediate. It is worth stressing that the halo 
mass interior to 5-Rd is almost five times that of the disc, yet its 
orientation does not follow that of the outer halo and massive torus, 
but instead follows that of the lighter disc, which precesses around 
the axis of the torus because of its angular momentum. Because 
the massive inner halo is oblate and aligned with the inner disc, 
it adds signific antly to the precess ion and curving of the LON in 
the outer disc. lBirmev et alj il998h also found that the inner halo 
quickly realigns with the disc, even when a large initial misalign- 
ment is imposed. 



5.2 Halo damping 

Figure fTsl shows that halo damping of the warp by the live halo is 
weak. Now for the convenience of description, we name the canon- 
ical run as "Run CL" (i.e. canonical run live halo). We introduce 
"Run CR" (i.e. canonical run rigid halo) with the almost same halo 
potential^ as in the canonical run except that the halo is rigid and 
stays exactly spherical. Panel (a) shows that nested rings of test 
particles on circular orbits in a rigid halo precess at constant incli- 
nation, dt, to the plane of the torus. Since rings precess at different 
rates, the apparent deviations from 15° at later times are caused by 
averaging rings with the same inclinations, yet different azimuths, 
over a finite radial range. (We have verified that the inclinations of 
individual particles with respect to the torus plane indeed remain 
close to 15°, as expected.) Panel (b) shows that self-gravity in a 
massive disc causes differential precession about the axis of the 
inner disc (Run CR). The quasi-periodic modulation of the inclina- 
tions of the outer rings in (b) is a manifestation of the extra torque 
from the inner disc. The initial decrease of the inclinations of most 
annuli to the torus plane (between f = and 600) has nothing to 

^ Note that the two potentials cannot be exactly the same, because the live 
halo is slightly oblate due to the gravity of the disc initially, and its oblate- 
ness evolves with time. 




Figure 12. The time evolution of the halo axes in the canonical run, (a) 
the tip angle and (b) ^lon both measured in the fixed initial frame. 
The solid curves, with decreasing thickness, represent the orientations of 
the inner halo (0 < r < 4_Rd), the central halo (4Rd < r < 8Rd) and the 
outer halo (8-Rd < f < 12_Rd) respectively, while the dashed curves show 
these angles for the inner disc. 



do with dynamical friction or halo damping since the halo is rigid 
in this case. The final panel, (c), shows the effect of using a live 
halo (the canonical Run CL). The precession due to the disc is aug- 
mented by the aligned flattening of the inner halo, which is partly 
responsible for the differences from panel (b), but a full explana- 
tion of all the details in this Figure is complicated; e.g., we would 
need to take into account the inclination of the halo at intermediate 
radii. The important result is that even after 2000 dynamical times, 
the mean inclination of the inner disc to the torus has not decreased 
by more than 3° ; i.e., damping from the live halo is very weak. 

O ur r esult differs from the conclu sions of Dubinski & Kuiikeg 
( 1995) and lNelson & Trema ine 1 1995) because the scenarios stud- 
ied are different. Their models have a very large initial misalign- 
ment between the halo and the inner disc and a significant pre- 
cession rate calculated according to the "modified tilt mode" in 
Soarke & Casertano ( 1988). Another difference is that the halos in 
their models were given no time to adjust to the strong and rapidly- 
precessing warp initially inserted in. 

In our simulations, the precession rate of the inner disc is 
very low, and the inner halo follows its motion closely, as shown 
in Figure [121 which is why dynamical friction between the two 
components is negligible. Only the outer halo is misaligned and 
damping comes mainly from coupling between the inner and outer 
halo. Such coupling is wea k for the isotropic halo models we used 
JNelson & Tremaindl 199511 . This is also consistent with the recent 
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Figure 13. (a) The time evolution of 8t, the tip angle with respect to the 
torus plane, of the inner disc (0 < R < 3-Rd , the black line) and 7 equal- 
width annuli spanning the radial range 3R^ < R < 7.9R^ (represented by 
the red, green, blue, cyan, magenta, yellow, and orange lines with increasing 
radius, respectively) for a rigid halo + test particle disc run. (b) As for (a), 
but for a massive, responsive disc in a rigid halo (Run CR). (c) As for (a), but 
for a massive disc in a responsive halo (the canonical Run CL). The slight 
deviations from 15° at later times in (a) are caused merely by averaging 
over a finite radial range. 



cosmological hydrodynamical simulations by iBailin et alJ |2] 
which found that the relative orientations of inner (r < O.lrvir 
and outer (r > O.lrvir) halo are uncorrelated. 



6 DISCUSSION 

6.1 The source of the external torque 

The torque in our model arises from a uniform, massive torus cen- 
tered on the galaxy and inclined to the disc, as originally adopted by 
Ijian g & BinnevI il99St) . The warp is driven by the quadrupole field 
of the torus, which varies as Mt/Rt, and therefore a less massive 
torus at a smaller radius has a similar effect, as we have verified in 
other experiments. However, a weaker perturbation would give rise 
to a milder warp that would take longer to form. 

A natural direct interpretation of the torus is an idealized rep- 
resentation of the stream of stars and dark matter from a disrupted 
orbiting satellite, if the mass is sufficiently large. There are numer- 
ous examples of accreted satellites, or companion galaxies orbiting 
in planes that are misaligned with that of the disc of the host galaxy. 
The Sgr Dwarf and the Magellanic Clouds are clear examples in the 
Milky Way, but their masses are generally on the low side unless 
their halos are substantial. 

On the other hand, the torque from the torus is similar to that 
of a flattene d outer halo, which may also have misaligned angular 
momentum I Ou inn & Binne V 19920. Conceptually, one could sep- 
arate the halo into two mass components: an outer halo flattened 
in a plane misaligned with the disc and having a uniform, quasi- 
spherical core, and the remainder - an inner halo with any arbi- 
trary radial mass profile, that may be strongly peaked to the centre 
of the galaxy. Regardless of how flattened the inner halo may be, 
it exerts no torque on the inner disc, since it is generally aligned 
with the disc plane (see i|5j- Iri this picture, the gravitational po- 
tential of the outer oblate halo alone might approximate the form: 
$oh(-R, 2) = Vo/2 ■ In (_Rc + R^ + z'^/ql), which arises from a 
flattened isothermal envelope with a quasi-spherical, uniform core 
of radius ~ _Rc , with a density scale set by vo . The ellipticity of the 
equipotential surfaces e# = 1 — gj. is approximately ep/3, where 
ep is the ellipticity of the halo density. It is easy to show the torque 
due to this outer halo exerted on a ring in the disc at radius r<Rc 
is (T(r)) ~ r^vo{<li ~ l)/(2^c)> which has the same radial de- 
pendence as the torque of a uniform torus. The effects discussed 
in i|4.4l and i|6.2l are dependent on the outward-increasing rate of 
retrograde precession of our model. The torque from our adopted 
torus is of the same order of magnitude as such an outer halo with 
a typical flattening in the range 0.6 < qp < 0.9 depending on R^. 

Cosmological simulations that include baryonic infall have 
not yet settled on a set of robust predictions for the shape and align- 
ment of the outer halo. A misaligned potential of the form just de- 
scribed is highly idealized, but has served as a useful theoretical 
exercise to understand warps. It has yielded some insight that may 
capture features of the torques that drive real warps. 



6.2 Persistence of warps 

The windup rate of a warp depends on the radial variation of the 
net precession rate. (The amount of the windup can be measured 
by cot i = AtR ■ dcop/dR, see lBinnev & Tremainel 19871 ') Forcing 
by our adopted torus causes differential precession of the opposite 
sign from that of the disc, making the total retrograde precession 
less differential, especially in the warp region (4_Rd < R< 10-Rd)- 
Thus the stiff inner disc and the two precessions work together to 
slow the windup of the warp. 

However, it is more interesting to study the evolution of a pre- 
warped massive disc without any external forcing. We therefore 
conducted two further experiments in which we removed the forc- 
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t = 700 




Figure 14. Removing the external forcing by the torus has little effect on 
the amplitude of the warp, but does affect the windup rate. The LON of Run 
CR (left) at t = 700 is slightly leading, but almost straight while that for 
Run CR-torus (right) at the same time is more wound up. Note that these 
runs, with a rigid halo, had the exactly same LON at i = 400. 



t = 700 




Figure 15. As Figure lT4l but for a live halo. The tip-LON plot from of Run 
CL (the canonical run) at t = 700, shown on the left, again has similar 
amplitude to that for Run CL-torus on the right but is less wound. 



ing torus after some evolution. We ran cases with both live (CL- 
torus) and rigid (CR-torus) halos, derived from the canonical runs 
CL (live halo) and CR (rigid halo). In both cases, we gradually re- 
moved all the torus particles in a linear fashion between t — 400 
and 600 after which, no torus particles remain. (Note that in all runs 
the torus mass was grown linearly to its maximum between i = 
and 200.) 

In these two runs, a warp was already well developed by t = 
400. Ramping down the torque from the torus allows us to study 
the unforced behaviour of the warped disc. In run CR-torus, the 
disc is affected by its self-gravity only (because the halo is rigidly 
spherical), but in run CL-torus, the effect of self-gravity of the live 
halo is also included. 

In Figure fT4l and Figure fTsl we compare the LON between runs 
in which the torus was left in place and the cases in which it was 
removed. Figure fT4l shows the two cases at i = 700 for the rigid 
halo, and Figure fTsI shows the corresponding two cases with the 
live halo. Note that both experiments in each pair are identical at 
t = 400. Two results are evident from this comparison, regardless 
of whether the halo is rigid or live: first, the amplitude of the warp 
is little changed by the removal of the forcing torus confirming that 
damping from the halo is weak; second, the LON in runs with the 
torus removed is more wound than in runs where the torus remains. 
This second point confirms that the torque from the torus is impor- 
tant to balance out the differential precession due to the stiff inner 
disc, and reduces the windup rate of the warp, but the extra winding 
is mild even ~ 3 Gyr after the torus starts to be taken away. Thus 
the warp can persist for a few Gyrs after the imposed quadrupolar 
perturbation is removed. 



6.3 Observations of the LON at large radius 

Our analysis has shown that the LON of warps in the inner region 
forms a leading spiral due to the potential of the disc and inner 
halo. The fact that the LONs of most warps form leading spirals 
over an extended radial range may also imply that the disc mass is 
significant in the central region of these galaxies (e.g. Figure|8}. 

The spirality of the LON at larger radii, on the other hand, 
contains information on the mass distribution in the outer halo that 
drives the warp. For example, if warps are formed by external per- 
turbations that exert torques of the form of our adopted torus, the 
LON should twist from leading to trailing at some large radius. The 
absence of a twist would imply a perturbation of a different form. 

We have looked for signs of a t wist in the largest sample of 
galaxies with tip-LON data to date iBriggs 1990). The transition 
radius _Rtw ~ 8-Rd for the torus we study ( i|4.4> . a value that has a 



mild dependence on the torus properties (Rt, 



m/Rtr°-\ 



If most of the large, late-type spirals in Briggs's sample have 
Holmberg radii in the range Ruo = R2e.5 ~ 4.4_Rd, we expect 
-Rtw ~ 2i?Ho and we therefore examine those having LON data 
beyond 2_Rho- 

Six galaxies (NGC 628, NGC 5055, NGC 2841, NGC 1058, 
NGC 3344 and M83) in this sample have LON data outside 2_Rho. 
We find that NGC 628 and NGC 5055 show signs of the LON spi- 
rality twist near 2i?Ho; NGC 2841 and NGC 1058 show weak signs 
of the spirality twist; the LON of M83 is quite straight near 27?ho, 
and might have spirality twist at an even larger radius; the LON of 
NGC 3344 is very noisy, but the spiral of the overall LON trend 
does not seem to twist. Thus perhaps four out of six galaxies reveal 
hints of spirality twists in their LONs. 

Thus currently available data do not provide convincing evi- 
dence for the spirali ty twist, mainly because the LON data beyond 
2i?Ho (Briggs 1990) are generally sparse and noisy. The outer mor- 
phology of the LON is worth more detailed investigation in deeper 
HI 21 cm observations of gas-rich galaxies, since it may contain 
valuable information on the shape and profile of the mass distribu- 
tion in the outer halo. 



6.4 The apparent thick disc due to warping 



The intrinsic thickness of the disc about its mean inclination is not 
greatly increased by warping. Comparison of the three panels in 
Figure fT6l reveals that the evolution of the disc thickness is little 
affected by an aligned torus, and a tilted torus causes only a slight 
increase for R > 4Rd. 

But a warped, luminous disc will contribute some light outside 
the thin, inner disc which, from some viewing directions (see Fig- 
ures^), may resemble a low-surface brightness thick disc. Sup- 
pose an idealized warp that starts from R25 at an angle a to the 
inner disc and has a straight LON. (We define a to be the kink an- 
gle between the flat inner disc and the steepest edge of the warp 
when seen in projection along the LON, i.e. it is not the warp 
angle 6 defined above.) The projected surface brightness profile 
along the minor axis of the edge-on disc will clearly be exponen- 
tial with vertical scale-height h^ — Rd sin a, and its brightness 
is roughly /(i?warp ~ R25)/ sin a ~ (25 + log^o sin q)/xb. The 
angle a ~ 45° for the warp in Figure [2^ and Figures [TTl Thus, if 
the warp starts in the luminous disc, it could possibly affect con- 
clusions about thick discs from photometry of edge-on galaxies. 
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Figure 17. The projected view of the warp in Figure|2ji when the viewpoint 
is rotated by 90° from that in Figurell^. 




radius 




radius 
Figure 16. Radial variation of disc thickness (rms z dispersion) relative to 
the mean inclination, at different times. The curves from bottom to top show 
values at t = 0, 100, 200, 300, and 400, respectively. The canonical simu- 
lation is shown in the top panel. The middle panel shows a comparison run 
with the same torus parameters as the canonical simulation except that the 
torus is in the disc plane, and the bottom panel shows another comparison 
run without a torus. 



7 CONCLUSIONS 

Now that the role of a live halo has been properly taken into ac- 
count, the idea that galaxy warps are manifestations of long-lived 
warp modes appears to have reached a dead end. But warps as re- 
sponses to external forcing remain viable, provided that suitable 
perturbations are frequent enough. Hierarchical galaxy formation 
scenarios guarantee late infall within gravitationally bound sub- 
groups, both diffuse and in lumps, and the infalling matter is un- 



likely to_shareJhe_an2ulaijnomento early disc and 

halo JOstriker & Binnevll989tlOuinn & Binnevll 992). 

The simulations presented here are motivated by this idea, but 
are limited to an in-depth study of one highly idealized case. We ex- 
amine how an initially flat particle-disc with random motion, em- 
bedded in a live halo, gradually acquires a warp as a result of a 
steady, applied external torque. Rather than striving for realism, we 
have endeavoured to understand how the warp develops and why 
the LON generally forms a loosely-wound, leading spiral. 

Despite having focused mainly on one case, we find a long- 
lived, large-amplitude warp that resembles those observed in many 
respects. The disc is flat in the inner part, and starts to warp at 
~ 5-Rd into an open leading spiral, which is consistent with the 
warping rules found by Briass ( 1990). 

The external torque causes the disc to precess in a retrograde 
sense. The warp develops because the precession rate of the inner 
disc, which is strongly cohesive due to its self-gravity and random 
motion, is slower than that of the outer disc. The growing misalign- 
ment between the inner and outer disc gives rise to a new torque 
acting on the outer disc from the massive inner disc, and the inner 
halo that is coupled to it. The torque from the interior mass is re- 
sponsible for the leading spiral of the line of nodes - our adopted 
external field alone would produce a trailing spiral. The fact that the 
LON of most warps forms a leading spiral over an extended radial 
range seems to imply a massive disc (e.g. Figure|8j. 

On the other hand, the leading spiral tells us little about the 
cause of the warp, but such information might be revealed by the 
behaviour of the warp at large radii. We show that the LON twists 
to trailing at very large radii from our experiments with a perturba- 
tion that alone would cause the warp to form a trailing spiral. Thus 
better observational data on the shape of the LON at very large 
radii may reveal the radial dependence of the torque that created 
the warp, and thereby provide information on the shape of the halo 
at large radii. 

Because the two separate components of the torque in our 
model cause differential precessions in opposite senses, the net ret- 
rograde precession of the warp is less differential than when only 
one is present. The moment we removed the external torque, the 
warp began to wind, although its amplitude did not decrease. Thus 
the quasi-steady warp we observe that lasts for more than a Hubble 
time is a consequence of steady external forcing. 

Even though the disc precesses due to the external torque, its 
motion is h ardly damped over many G yr, in contrast to the expec- 
tations from lNelson & Tremaineiri995i) . Damping is weak because 
the slow precession rate allows the inner halo to remain closely 
aligned with the disc, which therefore causes little drag. The very 
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weak damping seems to be caused more by the relative precession 
of the inner and outer halo. 

A fixed outer torus is clearly unrealistic and external forcing is 
likely to be strongly time-dependent in reality. We cannot say much 
about the rate at which misalignments in real halos will settle, but 
we have shown that the warp survives for a few gigayears after the 
torus is removed. Time dependence may anyway be a side issue if 
the h alo axis is constantly slewing, as argued bv lOuinn & BinnevI 
il992) . Warps formed this way can be repeatedly regenerated when 
a new infall event happens. Since cosmic infall and mergers are 
more likely to happen in a denser environment, warps can be in- 
duced more frequently in such an environment, wh ich is consistent 
with the warp statistics in Garcia-Ruiz et al. ( 2002). 

Since warps are ubiquitous and a gravitational phenomenon 
closely related to hierarchical galaxy formation, their properties 
may possibly be able to constrain small-scale cosmic structure. 
Clearly further experimentation with various external perturba- 
tions, halo profiles and masses, disc properties (random motions, 
barredness), etc., is required to discover what might be inferred 
from the existence of warps in galaxies. 
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ring accreting torus 

Figure Al. The calculation of the precession rate of a disc annulus, at radius r and inclined relative to the torus symmetry plane with the angle i, under the 
torque of the accreting torus. The precession is retrograde about the symmetry axis of the torus. 

APPENDIX A: THE PRECESSION RATE DUE TO THE TORUS 



Suppose the disc is composed of many rigid rings (or annuli). Here we calculate the retrograde precession rate of one such ring, at radius r 
and inclined relative to the torus symmetry plane with the angle i, due to the potential of the accreting torus (see Figure IaTI . 
The potential of the accreting torus (its thickness ignored) can be obtained as the following, by solving Laplace's equation. 



$(r, 



' Rt 



^ ~ I (:^)'^'^'°'^) + I (;J) V4(COS0) + (rVi?? 



(Al) 



where 6 is the conventional polar angle in the polar coordinate system based on the torus; Rt and Alt are the radius and mass of the torus, 
respectively. 
We define 



r(r) 



1 GMt 



and Q{r) 



3 GMt / r \ 4 



2 Rt \ Rt ) 8 -ft 

for convenience, so to the order of r'^ /Rt we have 
GMt 



t \Rt) 



$(r, e) = 



Rt 



+ T(r)P2(cos6l)- Q{r)Pi{cos0). 



The only force component relevant to the calculation of torque on the annulus is Fe : 

r oO r dcosO r dcosO 

Since 



(A2) 



(A3) 



(A4) 



P.(.) = -(3x -1); -^^ 



3x: 



Pa{x) = i(35x* - ZQx^ + 3); and ^^^ = i(35x^ - 15x), 

we find 
Fb = 3^^ cose sin( 



iQW 



2 r 



(15 cos ^ 



35 cos 9) sin ( 



(A5) 



The torque exerted on a point mass in the rigid ring due to the torus depends on the azimuth of the point mass; the maximum value of 
the torque occurs at azimuth = or vr where 9 — -k 12 — i. 



|r X F| =rF6 ^3T{r)cos9 shi9 + -Q{r){15 - 35cos'^ 9)cos9 sm9 

1 2 

3r(r)sini cosiH — (3(r)(15 — 35sin i)sini cosi. 



{r) 



The azimuthally average torque on the rigid ring is therefore 

p2TT 

_ Jo 



: COS^ (j) I 



13 1 

Tmax = -r(r)sinj cos j + —Q(r) (15 — 35 sin i)sini cosi. 



27r 2 2 

We can obtain the total retrograde precession rate due to the torus (up to the fourth order) 



(A6) 



(A7) 



(r) 3 T(r) . 1 g(r)(15 - 35sin2 i) 

UJp = -r-. : = TJ- COSI + — cos t 

Lsim 2r -Vc 4 r ■ Vc 



3 GMtr 



4 7??v; 
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C O' 




Figure Bl. The schematic diagram showing how and 4"lon are calculated. 



,(2) 



i(^)^15-35sin^z) 



For qualitative purposes, we generally use the following simpler second-order approximation (ujp ), unless the radius of the annulus is 



(2)^, 



(A8) 



large. 



,(2) 



3 GMtr 



(A9) 



" P 4 i??K 

Note that same results can be obtained for the nodal precession rate of a test particle's circular orbit (as opposed to a rigid ring), 
the only difference here is that th e average is taken as the time average over one orbit (the orbit-average procedure can be found in 
ISteiman-Cameron & Diirisenll984. Appendix). 



APPENDIX B: THE DERIVATION OF TIP ANGLE OF RINGS AND AZIMUTHAL ANGLE OF THEIR LON 

We first calculate the tip angle O relative to the original z-axis. We know AB = 20B sin y . On the other hand we also have 

AB ^20'B sin— ^ =20Bsmi sin— ^ (Bl) 

here A$p = J LUp{t) dt is angle that the annulus has precessed. So obviously the equation to determine Q is 

sin — = sin i sin — — -£■ (B2) 

If we define the azimuthal angle of the LON of a ring as the counter-clockwise angle from -\-x axis, "3?lon = /.ADC — n. It can be 
shown, with some trigonometry on triangles in Fieure lBTI that 



lADC = cos"^ 



1 — sm I sm -^ 



(B3) 
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So 



"l>LON = lADC ~ -K — COS ^ 



1 - sin2 i sin2 ^ 



$LON 



For small i, we have 



(B4) 



(B5) 
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